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The present paper extends to clusters of galaxies the study of Del Popolo (2012), 
concerning how the baryon-dark matter (DM) interplay shapes the density profile of 
dwarf galaxies. Cluster density profiles are determined taking into account dynamical 
friction, random and ordered angular momentum and the response of dark matter 
halos to condensation of baryons. We find that halos containing only DM are char- 
acterized by Einasto's profiles, and that the profile flattens with increasing content 
of baryons, and increasing values of random angular momentum. The analytical re- 
sults obtained in the first part of the paper were applied to well studied clusters whose 
inner profiles have slopes flatter than NFW predictions (A611, A383) or are character- 
ized by profiles in agreement with the NFW model (MACS J1423.8+2404, RXJ1133). 
By using independently-measured baryonic fraction, a typical spin parameter value 
A ~ 0.03, and adjusting the random angular momentum, we re-obtain the mass and 
density profiles of the quoted clusters. Finally, we show that the baryonic mass inside 
~ 10 kpc, Mb^n is correlated with the total mass of the clusters, as Mb,i n oc M®qq. 
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1 INTRODUCTION 

The ACDM model is remarkably successful in fitting a wide 
range of data of large-scale structures, starting from the 
anisotropy and polarization spectrum of cosmic microwave 
background radiation (Spergel et al. 2003; Komatsu et al. 
2011), passing through the high-redshift studies of the Lya 
forest, and going on with the Hubble diagram of Type la Su- 
pernovae (Kowalski et al 2008), and the matter power spec- 
trum with its Baryonic Acoustic Oscillation features (Perci- 
val et al. 2010). 

In spite of the success of ACDM predictions on 
cosmological scales, the model has shown difficulties in 
giving correct predictions on galaxy scales, especially for 
what concerns the central density expected in dark matter 
haloes. Dwarfs and low surface brightness (LSB) galaxies, 
objects that are dark matter dominated, are characterized 
by kinematics, and rotation curves, incompatible with 
haloes predicted by ACDM, as obtained in dissipationless 
N-body simulations (Navarro et al. 1996, 1997, hereafter 
NFW; Moore et al. 1998; Jing & Suto 2000, Klypin et al. 
2001; Hayashi et al. 2004; Power et al. 2003; Navarro et al. 
2004, 2010). 



E-mail: antonino.dclpopolo@unibg.it 



While dissipationless N-body simulations predict cuspy 
profiled, observations show that the inner part of density 
profiles is characterized by a core- like structure (e.g., Flores 
& Primak 1994; Moore 1994; Kravtsov et al. 1998; de Blok & 
Bosma 2002; de Blok, Bosma & McGaugh 2003; Gentile et 
al. 2004, 2006; Blaise-Ouelette et al. 2004, Spano et al. 2008, 
Kuzio de Naray et al. 2008, 2009, and Oh et al. 2010). The 
quoted discrepancy between observations and simulations is 
known as the Cusp- Core problerrJ3. 

Another important test of the results of numerical 
simulations is to check their predictions also on cluster 
scales. Constraints on DM density profiles are more 
straightforward to obtain in clusters of galaxies than in 
galaxies0 since clusters have several measurable properties 
that can be understood and interpreted in a simpler fashion 



1 The inner density profile predicted by simulations is character- 
ized by p oc r _1 for the NFW model, p oc r~ 1,5 in the Moore et 
al. (1998) model, and by the Einasto profile, shallowing towards 
the center of the halo, in more recent simulations (e.g., Navarro 
et al. 2010). Stadcl et al. (2009) found a ~ 0.8 at 120 pc. 

2 For the sake of precision, we have to report that some obser- 
vations found density profiles compatible with both cuspy and 
cored profiles (van den Bosch et al. 2000; Swaters et al. 2003a, b; 
Spekkens et al. 2005; Simon et al. 2005; de Blok et al. 2008). 

3 Even if a drawback in using clusters is that the mass distribu- 
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than galaxies rotation curves. For example, X-ray emission 
from the intracluster plasma at radii of the order of 10% 
of the virial radius (r v ir) can easily observed and the gas 
temperature profile can be measured with observations 
(e.g., with Chandra or XMM). 

Nevertheless many studies have used X-ray observa- 
tions of the hot intracluster medium, under the assumption 
of hydrostatic equilibrium, X-ray data alone have difficulties 
in constraining the mass distribution, especially in the cen- 
tral regions. This because even relaxed clusters tend to have 
"cooling flows". In these clusters X-ray emission is often 
disturbed and the assumption of hydrostatic equilibrium is 
questionable (see Arabadjis, Bautz & Arabadjis 2004). 

X-ray temperature measurements are carried out 
from 500 kpc (Bradac et al. 2008) to ~ 50 kpc, and the 
determination of temperature at smaller radii are limited by 
instrumental resolution or substructure (Schmidt & Allen 
2007). X-ray analyses have obtained wide ranging values 
of the inner slope, a, of the density profile with a ranging 
from ~ 0.6 (Ettori ct al. 2002) through ~1.2 (Lewis, Buote 
& Stocke 2003) to ~1.9 (Arabadjis, Bautz & Garmire 
2002), while Chandra and XMM-Newton results suggest 
good agreement with CDM predictions (see Schmidt & 
Allen 2007, and references therein). The value of a obtained 
by using X-ray observations have the drawback that it is 
complicated to take account of the stellar mass contained in 
the brightest cluster galaxy (BCG), located in the cluster 
center. The central mass, constituted by stars and gas, even 
if it is small, compared to the total mass of the system, 
usually dominates the mass at small radii (~ 10 kpc), with 
strong implications on the shape of the inner DM profile. 

The dark matter distribution in clusters can be 
also studied through gravitational lensing. Weak lensing 
of background galaxies is used to reconstruct the mass 
distribution in the outer parts of clusters (Mellier 1999). 
This technique is based on averaging noisy signal coming 
from many background galaxies. The resolution that can be 
achieved is able to constrain profiles inside ~100 kpc. An 
example of its application was given by Dahle, Hannestad 
& Sommer-Larsen (2003), who found that the average 
profile of 6 clusters agrees with the NFW profile at radii 
r ^ 0.1r„i r , even though with large uncertainties. 

In the central parts of the cluster, lensing effects 
become non-linear and one can use the strong lensing 
technique, in order to constrain the mass distribution. This 
technique has a typical sensitivity to the projected mass 
distribution inside ~100-200 kpc, with limits at ~10-20 kpc 
(Gavazzi 2005; Limousin et al. 2008). Typical structures 
observed in the strong lensing regime are radial arcs, 
located in positions corresponding to the local derivative 
of the cluster mass density profile, and tangential arcs 
whose position is determined by the projected mass density 
interior to the arc. 



tion in clusters is usually affected by the growth of the central 
galaxy. 



Example of application of the lensing method are Smith 
et al. (2001), who studied the tangential and radial arcs 
of A383, finding an inner slope of the total mass density 
profile a > 1 at a radius ~ l%r V i r . Sand et al. (2004), 
and Newnman et al. (2011), after removing the baryonic 
component, with the aid of stellar kinematics, found that 
the dark matter only inner slope is flatter than 1. Kneib et 
al. (2003) studied CI 0024+1654 by means of strong and 
weak lensing finding that a NFW profile fits the profile from 
0.1r„i r to several values of r v w (but Tyson, Kochansky 
& Dell' Antonio 1998 found a = 0.57 ± 0.024). Gavazzi 
et al. (2003), and Gavazzi (2005) studied MS2137.3-2353 
concluding that the precise value of the slope depends on 
the mass-to-light ratio of the BCG (but Sand et al. 2002, 
found a cored profile with a ~ 0.35). 

As it is clear from the previous examples, strong 
gravitational lensing has given controversial results, and 
in several cases a values are much smaller than ACDM 
predictions. Moreover, similarly to X-ray analyses, lensing 
alone cannot disentangle the DM and baryonic components. 
Better constraints on the central part of the density profiles 
can be obtained through stellar kinematics of the central 
galaxy (~l-200 kpc region), and even better results are 
obtained combining all the previous methods (Miralda- 
Escude 1995; Kneib et al. 2003, (weak+strong lensing); 
Bradac et al. 2005 (weak+strong lensing); Mahdavi et al. 
2007 (X-ray+weak lensing)). 

As previously reported, distinct researchers, using dif- 
ferent techniques, found results going from good agreement 
with the ACDM model to disagreement with it, and what is 
worse is that this often happened for the very same clusters 
(Smith et al. 2005; Gavazzi et al. 2005; Zappacosta et al. 
2006; Schmidt & Allen 2007; Umetsu & Broadhurst 2008; 
Bradac et al. 2008; Limousin et al. 2008). 

The quoted discrepancies in observational results are 
probably connected to several factors: a) use of observa- 
tional techniques with different/limited dynamic range in 
radius; b) different definition of the slope, which sometime 
refers to the DM and sometime to the total mass. We 
should also not forget the scatter in the profile from cluster 
to cluster, also observed in the case of dwarfs by Simon 
et al. (2005), which contributes to the quoted discrepancy 
between observational results; c) not taking into account 
the stellar mass of the BCG, which is correct for the outer 
part of the cluster but wrong for the inner part, where stars 
dominate the density profile. Concerning this last issue, if 
the baryon-DM interaction changes the density profile, the 
usual parameterization of the density profile with NFW or 
generalized Navarro-Frenk- White (gNFW) models, may be 
inappropriate (Sand et al. 2008). Another important factor 
that can produce errors is originated if the clusters shape is 
not taken correctly into account (Gavazzi 2005). 

A series of papers by Sand (Sand et al. 2002, 2004, 
2008, hereafter S02, S04, and S08), further fed the debate 
on the shape of DM halos of clusters, and the inner slope of 
the density profiles. In the quoted papers, Sand combined 
gravitational lensing with the velocity dispersion profile of 
the BCG. In this way the author was able to separate the 
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contribution to the halo coming from the DM from that 
coming from the stellar mass of the BCG. 

S02 studied MS 2137-23 finding a flat inner slope with 
a < 0.9 at 99% CL, and in their best fit parameters, a was 
0.35. 

S04 studied a sample of 6 clusters (MS 2137-23; A383; 
A963; RXJ1133; MACS 1206; A1201), three having just 
tangential arcs, and three also having radial arcs. Com- 
bining lensing and velocity dispersion profile of the central 
galaxy in each case, they found a mean DM distribution 
inconsistent with the NFW value, a = 1 at > 99% CL. The 
system of clusters containing radial arcs gives a value of 
a = 0.52±g;g| (68% CL), and a < 0.57 at 99% CL, for the 
case of the tangential arc sample. If their results are correct, 
they disagree with ACDM model on small scales, unless 
some physical mechanism not taken into account could get 
the ACDM out of the trouble (e.g., a form of interplay 
between DM and baryons). Bartelmann & Meneghetti 
(2004) and Meneghetti et al. (2007) concluded that the 
flat slopes obtained by Sand were due to an oversimplified 
assumptions in the analysis, such as negligible ellipticity. 

An improved analysis of A383 and MS2137-23, 
presented by S08, found a flat value for the slope of 



A383 (a = 0.45 



+0.2 n 



Newmann et al. (2009) (hereafter 



N09), presented a detailed analysis of DM and baryonic 
distribution in A611, combining weak lensing, strong 
lensing and stellar velocity dispersion for the BCG, finding 
a slope a < 0.3 (68% CL) and similarly, Newman et 
al. (2011) (hereafter Nil) found a < 1 (95% CL) for 
A383. They showed that degeneracies in constraining 
the DM profile can be broken only simultaneously using 
the three techniques. The quoted results show that at 
least some clusters of galaxies have inner density profile 
slopes shallower that those obtained in N-body simulations. 

These results are in agreement with studies on dwarf 
galaxies, showing that not all dwarfs have core-like rotation 
curves and that their density profiles are compatible with 
both cuspy and core profiles (Hayashi et al. 2004; van den 
Bosch et al. 2000; Swaters et al. 2003a,b; Simon et al. 2005, 
Spekkens, Giovanelli & Haynes 2005; de Blok et al. 2008). 
Very interesting is the case studied by Simon et al. (2003, 
2005) who found that in a sample of dwarf galaxies, namely 
NGC 2976, NGC 4605, NGC 5949, NGC 5963, and NGC 
6689, the inner slopes goe from very flat to cuspy, more 
precisely a = (0.01; 0.78; 0.88; 1.20; 0.79), respectively, with 
a mean slope a ~ 0.73, and a dispersion of 0.44. 

de Blok et al. (2008), using a sample from the HI 
nearby galaxy survey (THINGS) found that the best fit to 
rotation curves depends on their mass: for galaxies having 
Mb < —19 the NFW profile or the pseudo-isothermal 
profile statistically fitted equally well, while this what not 
the case galaxies with Mb > —19 the core-dominated 
pseudo-isothermal model fitted significantly better than the 
NFW model. 

In several previous papers, we studied the role of 
baryons on the previously quoted problem (e.g., DP09, 
DP12). 

In the present paper, we extend the analysis of DP12 



from dwarf galaxies to clusters of galaxies, studying the ef- 
fects of baryons on the inner DM density profile. Similarly 
to DP12, the present paper has two aims. The first is to 
study the role of baryons and random angular momentum 
in changing clusters of galaxies density profiles. The second 
is to apply the results of this study to observed clusters, to 
understand if the model can re-obtain the density profile of 
some well studied clusters with cuspy, cored and intermedi- 
ate density profiles. 

This paper is organized as follows. Section 2 summa- 
rizes the model used. In Section 3.1, we find how changes 
in baryonic fraction and angular momentum, influence the 
density profiles In Section 3.2, we compare the results of 
Section 3.1 with four different clusters of galaxies. Section 4 
is devoted to discussions and Section 5 to conclusions. 



2 SUMMARY OF THE MODEL 

The density profiles of clusters are obtained using the model 
discussed in DP09, and DP12. For convenience of the reader, 
the model is summarized in the following of this section. 

DP09 is an improvement of previous spherical infall 
models (SIM) presented in literature (Gunn & Gott 1972; 
Fillmore & Goldreich 1984; Bertschinger 1985; Hoffman & 
Shaham 1985; Ryden & Gunn 1987; Avila-Reese, Firmani 
& Hernandez 1998; Subramanian et al. 2000; Ascasibar, 
Yepes & Gottleber 2004; Williams, Babul & Dalcanton 
2004). 

Differently from previous SIMs, DP09 includes in the 
model the joint effects of ordered and random angular mo- 
mentum, dynamical friction, adiabatic contraction of dark 
matter. Previous SIMs usually took account one effect at a 
time: random angular momentum (e.g., Williams, Babul & 
Dalcanton 2004), dynamical friction of stellar/DM clumps 
against the background halo (e.g., El-Zant et al. 2001; 
Romano-Diaz et al. 2008), or just adiabatic contraction 
(e.g., Blumenthal et al. 1986; Gnedin et al. 2004; Klypin, 
Zhao, and Somerville 2002; Gustafsson et al. 2006). 

The main features of the SIM are described by Gunn & 
Gott (1972). The quoted model assumes that a protostruc- 
ture is divided into mass shells, each one expanding with 
the Hubble flow from an initial comoving radius Xi to a 
maximum one x m (usually named turn-around radius xt a ), 
and then collapse. Non-linear processes convert the kinetic 
energy of collapse into random motions, giving rise to a 
"virialized" structure. The density profile at turn-around, 
pta(x m ), is obtained following the dynamics of collapse 
of the shells, assuming that mass is conserved and that 
each shell is kept at its turn-around radius (Peebles 1980; 
Hoffman & Shaham 1985; White & Zaritsky 1992). 

In order to go on from turn-around to the final collapse, 
and to obtain the final density profile, one usually assumes 
that the potential well near the center varies adiabatically 
(Gunn 1977, Fillmore & Goldreich 1984). The final density 
is given by 



p{x) = 



Pta(Xm) 



1 + 



din/ 



d In x,- 



(1) 
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where / = x/x m is the collapse factor (see Eq. A18, DP09). 
During expansion, the tidal fields originated by the large- 
scale structure exert a torque on the proto-structure im- 
parting angular momentum on the proto-halo (Hoyle 1953; 
Peebles 1969; White 1984; Ryden 1988; Eisenstein & Loeb 
1995; Catelan & Theuns 1996). After that the proto-strucure 
decouples from the background, turns-around and collapses, 
tidal torque is noteworthy reduced, since the lever arm is re- 
duced drastically by collapse (Schaefer 2009). Angular mo- 
mentum originated by tidal torques is calculated obtaining 
the rms torque, r(r), on a mass shell and then calculat- 
ing the total specific angular momentum, h(r,v), acquired 
during expansion by integrating the torque over time (Ry- 
den 1988, Eq. 35). As shown in DP09, the values obtained 
for angular momentum are in agreement with others studies 
(e.g. Catelan & Theuns 1996). As in DP12, total angular 
momentum is expressed in terms of the dimensionless spin 
parameter 



A 



L\E\ 



1/2 



GM 5 / 2 ' (2) 
where L is the angular momentum, summed over shells, 
and E is the binding energy of the halo. According to Pad- 
mannabhan (1993), A can be interpreted as the ratio of the 
angular velocity, ui of the system to that, u) S u P of the system 
providing the rotational support, obtaining 

A = — = L (*\ 

oj sup 2G 1 / 2 M*/ 2 R 1 / 2 ' y ' 

In a system constituted by DM and baryons as ours, we have 
a spin parameter for DM and another for gas 



A 



gas(DM) 



^gas(DM) 



M, 



gas(DM) 



[2G{Mgas + M DM )r v 



1/2, 



(4) 



where L 



gas(DM) 



is the angular momentum of gas(DM), and 



gas(DM) 



the gas(DM) mass inside the virial radius r v i 



The ratio of X gas /^DM is fixed according to Gottleber & 
Yepes (2007) (1.23 for haloes with M vir > 5 x 10 14 /i _1 M Q 
(see their figure 5)). The distribution of the A parameter 
is well described by a lognormal distribution (e.g. Vivitska 
et al. 2002). Gottleber & Yepes (2007) studied the spin pa- 
rameter of 10,000 clusters extracted from the Mare Nostrum 
Universe SPH simulation, finding a log-normal distribution 
with best fit parameters given by A with A = 0.0351±0.0016, 
a x = 0.6470 ± 0.0067 for the DM distribution and A = 
0.0462 ± 0.0012, cr\ = 0.6086 ± 0.0030 for the gas distribu- 
tion, with A max = 0.0231, 0.0319 for DM and gas distribu- 
tions (see also Bett et al. 2007; Sharma & Steinmetz 2005). 
A spin parameter of A ~ 0.05 is characteristic of objects 
having negligible rotational support, and little systematic 
rotation (e.g., large elliptical galaxies and clusters of galax- 
ies). Spiral and LSB galaxies are rotationally supported and 
in the case of LSBs 0.06 < A < 0.21 are LSBs (Jimenez et 
al. 1998; Boissier et al. 2003). 

Apart from this source of angular momentum, origi- 
nated by bulk motions and tidal fields, a random angular 
momentum, j, is generated by random velocities (see Ryden 
& Gunn 1987). In several of the papers which took account 
of angular momentum in SIM (e.g., Nusser 2001, Hiotelis 
2002; Williams et al. 2004) only random angular momentum 
was taken into consideration. The approach usually followed 
was that of assigning a specific angular momentum at turn 



around (e.g., Nusser 2001; Hiotelis 2002; Ascasibar, Yepes 
& Gottleber 2004), having typical value: 



j = j» oc V GMr„ 



(5) 



The specific angular momentum, j, may be also expressed 
through the ratio e = {j^^j , where r mi „ and r max , 
are the pericentric and apocentric radii. Avila-Reese et al. 
(1998) used eo as a free parameter to take into consideration 
processes related to mergers and tidal forces that could pro- 
duce tangential perturbations. In their paper, they showed 
that the detailed description of these processes is largely 
erased by the virialization process, remaining only through 
the value of eo, which they fixed to eo = 0.3. The value 
e ~ 0.2 gives density profiles very close to the NFW profile 
(Avila-Reese et al. 1998, 1999). The previous procedure is 
based on results of N-body simulations of CDM halo col- 
lapse, giving constant < J.'"'" >~ 0.2 ratios of dark matter 
particles in virialized haloes. An improvement to the previ- 
ous model is due to Ascasibar, Yepes & Gottleber (2004), 
who found that particle orbits are slightly more radial as 
we move out to the current turn around radius, rt a , and 
that there is a dependence on the dynamical state: major 
mergers are well described by constant eccentricity up to 
the virial radius, relaxed systems are more consistent with 
a power-law profile, while minor mergers are in the middle. 
A least-square fit to the relaxed population yields: 



e(rmax) ^ 0-8(rmax/rta) 



(6) 



for r max < 0.1r ta . 

In the present paper, random angular momentum is 
taken into account through Avila-Reese et al. (1998) method 
with the correction of Ascasibar, Yepes & Gottleber (2004) . 

Dynamical friction was taken into account as described 
in Antonuccio-Delogu & Colafrancesco (1994) (see also Ap- 
pendix D of DP09). Its effects on structure formation were 
obtained by introducing the dynamical friction force in the 
equation of motion (Eq. A14 in DP09). We also took into 
account adiabatic contraction (AC) of dark matter halos in 
response to the condensation of baryons in their centers, and 
leading to a steepening of the dark matter density slope. 

We summarize in the following the AC model taking 
also account of angular momentum exchange. 

The present standard form of the AC model was intro- 
duced and tested numerically by Blumenthal et al. (1986) 
(see also Ryden & Gunn 1987; Barnes 1987; Oh 1990). Let's 
consider a spherically symmetric protostructure that con- 
sists of a fraction of baryons, Fb = Mb/Msoo << 1, and a 
fraction 1— Fb of dark matter particles constituting the halo. 
Here, M500 is the mass contained within a radius where the 
enclosed density is 500 times the critical density of the Uni- 
verse, and Mb = M* + M gas , where M* is the stellar mass 
and Mgas is that of gas. Baryons and halo particles are as- 
sumed to be well mixed initially (i.e., the ratio of their den- 
sities is Fb through the protostructure). As the baryons dis- 
sipatively cool and fall into a final mass distribution Mb{r), 
a dark matter particle initially at radius r\ will move in to 
radius r < ri, characterized by 



r [M b (r) + M dm (r)] = nMi(n) 



(7) 



where Mi(ri) is the initial total mass distribution, Mb(r) 
is the final mass distribution of dissipational baryons and 
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Mdm is the final distribution of dissipationless halo particles. 
Assuming that the orbits of the halo particles do not cross, 
one obtains 



M dm {r) = (l-F b )Mi(n) 



(8) 



Eqs. (0, ® can be iteratively solved to calculate the fi- 
nal radial distribution of the halo particles once Mi(ri) and 
Mb(r) are given. A usual assumption is that initially baryons 
had the same density profile as the dark matter (Mo et al. 
1998; Cardone & Sereno 2005; Treu & Koopmans 2002; Kee- 
ton 2001), and the final baryons distribution is assumed to 
be a Hernquist configuration (Rix et al. 1997; Keeton 2001; 
Treu & Koopmans 2002). In the present paper, the Blumen- 
thal's model was improved following Gnedin et al. (2004) 
who proposed a simple modification which describes numer- 
ical results more accurately. They proposed a modified adi- 
abatic contraction model based on conservation of the prod- 
uct of the current radius and the mass enclosed within the 
orbit-averaged radius: 



M(r)r = const. 

where the orbit-averaged radius is 



2_ 



dr 



(9) 



(10) 



where T r is the radial period. 

The previous, AC model, assumes no angular momen- 
tum exchange between different components (e.g., baryons 
and DM). In the early phase of collapse, the baryons den- 
sity is an order of magnitude smaller than that of DM, and 
the exchange of angular momentum has little impact on the 
dark matter. During the collapse baryons infall into the cen- 
ter of the proto-structure deepening the gravitational po- 
tential well. In the late stage of collapse the baryon density 
becomes large, non-axisymmetric component may develop 
due to the excitation of spiral waves and/or bar-like modes. 
The resulting increase of DM density acts as a sort of cou- 
pling process of the baryons and the DM (Klypin et al. 2001; 
Klypin, Zhao, and Somerville 2002). At later stages of col- 
lapse, when baryon density increases, the effect is potentially 
quite large resulting in a decrease of the dark matter density 
by a factor of ten (Klypin, Zhao, & Somerville 2002). 

In order to describe quantitatively the quoted process, 
let's consider a spherical shell of dark matter with density 
Pdm, radius r, thickness dr, and specific angular momentum 



j = rV c = y/G [M b (r) + M dm {r)]: 
one can get an implicit equation for the final radius 77 
AAM 



Jf = J 



1 + 



47rpd m r 3 



A = 1 + vf 

V c dr 



AM 



M h ( - M b . 



(11) 

(12) 

(13) 
(14) 



Klypin, Zhao, and Somerville (2002), where Mb,f is the final 
baryons mass and M — Mdm + Mh is the total mass inside 
a radius r. Eq. (|12|) is solved numerically. The solution also 
gives the mass inside a final radius 77. Eq. (|12[) has the 
same structure as Eq. (|7|), the only difference is the term on 
the right-hand-side, which is the correction due to angular 
momentum deposition. 



For what concerns baryons, we discussed how they were 
introduced in appendix E of DP09, and how their initial 
mass distribution and the final distribution were fixed. The 
quantity of baryons used in the calculation of Sect. 3.1, was 
fixed using Giodini et al. (2009) results. As shown from Mc- 
Gaugh et al. (2010) (see their Fig. 2), the fraction of baryons 
detected in all forms deviates monotonically from the cosmic 
baryon fraction as a function of mass. On the largest scales 
of clusters, most of the expected baryons are detected, while 
in the smallest dwarf galaxies, fewer than 1% are detected. 
The detected baryon fraction, fd, is the ratio of the baryon 
fraction, Ft = Mj/Msoo, and the universal baryon fraction, 
fb, is 



U = (M b /M 500 )/f b = F b /f b 



(15) 



where f b = 0.17 ± 0.01 (Komatsu et al. 2009). Note that 
in this paper the masses, M200, and M v i r , are converted 
to Msoo using the method used in White (2001), Hu & 
Kravtsov (2003), and Lukic et al. (2009). 

In the four clusters studied in Sect. 3.2, we calcu- 
lated the baryon fraction using Schmidt & Allen (2007) and 
checked that the result was in agreement with Giodini et 
al. (2009) results. More in detail, in this last case, we cal- 
culated the baryon content subtracting the total mass of 
clusters, obtained from the data in table 3 of Schmidt & 
Allen (2007), from the DM mass obtained from the data in 
their table 4, namely 



M b = M 



DM+b 



M L 



8007T 3 

7, Pc\ r vir,DM+ b 



r dm) 



(16) 



where p c is the mean background (critical) density, and 
r V ir = cr s , where the concentration parameter c, and the 
scale factor r a are given in Schmidt & Allen (2007) for 34 
massive, dynamically relaxed clusters of galaxies. Note, 
that Eq. (|16|l comes from the fact that for a NFW profile 

Pvi{r — Tvirj — 3 pcr v i r - 

In order to have a more clear view of how the previous 
quoted processes works to build up the cluster, we sum- 
marize the structure formation starting from high redshifts. 
Initially the proto-structure is in the linear phase, it ex- 
pands, reaches a maximum of expansion and then collapses. 
Baryons trapped inside the potential wells of DM halos are 
subject to radiative dissipation processes which give rise to 
clumps and self-gravitating clouds before it collapses to the 
halo center and condenses into stars and galaxies. The stage 
of baryons cooling and stars formation happens as described 
in Ryden (1988) (Sect. 4: "baryionic dissipation"). In the in- 
fall, the baryons compress the DM halos (AC), producing a 
steepening of the DM density profile. When clumps reaches 
the central high density regions they experience a dynami- 
cal friction force from the less massive DM particles as they 
move through the halo. Dynamical friction acts as an an- 
gular momentum engine (Tonini, Lapi & Salucci 2006), and 
energy and angular momentum are transferred to DM (El- 
Zant et al. 2001, 2004), increasing its random motion, and 
giving rise to a motion of DM particles outwardfl More- 

4 Another model of baryon-DM interaction, and exchange of 
angular momentum between the baryons and DM was previously 
discussed in this section (Klypin et al. 2001; Klypin, Zhao, & 
Somerville 2002). 
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over, ordered angular momentum mainly acquired in the 
expansion phase, through tidal torques, gives rise to nonra- 
dial motions in the collapse phase that amplifies the effects 
of the previous mechanism. Then the joint effect of angu- 
lar momenta (ordered and random) and dynamical friction, 
transporting it from baryons to DM, overcomes that of the 
AC and the profile starts to flatten. 

As previosly described, in order to follow the structure 
formation, the protostructure is divided into mass shells, 
containing DM and baryons, which after radiative dissipa- 
tion forms clumps, stars, and galaxies. In order to form 
the structure the evolution in various phases is followed by 
means of the SIM. After turn-around shells move inward, 
baryons (namely, baryonic clumps) exchange angular mo- 
mentum with DM and accretes to the center, forming the 
central galaxy (BGC) (see also Lackner & Ostriker 2010). 

As already shown in DP09, and DP12, the density pro- 
files, obtained with the model summarized in this section, 
and fully described in DP09, DP12, produces profiles in 
agreement with other studies (e.g., El-Zant et al. 2001, 2004; 
Mashchenko et al. 2005, 2006, 2007; Romano-Diaz et al. 
2008, and the more recent paper of Governato et al. 2010). 



3 RESULTS AND DISCUSSION 

3.1 Changes of density profiles with baryons 
fraction and angular momentum 

As discussed in the introduction, dwarf galaxies and LSBs 
have rotation curves incompatible with the haloes that A 
CDM model predicts. This problem, known as Cusp-Core 
problem, pointed out by Flores & Primack (1994), and 
Moore (1994) has not seen till now a clear-cut solution. 

Inner density profiles in clusters of galaxies have also 
been at the center of a similar debate to that of dwarfs, 
with claims of value of a going from a = 0.35 (S02) to ~1.9 
(Arabadjis, Bautz & Garmire 2002). 

To this aim, in the following of the paper, we use DP09 
to generate density profiles in the range 1O 14 -1O 15 M , then 
focalizing our attention on changes of inner density profiles 
produced by different values of baryon fraction and random 
angular momentum. 

In Sect. 3.1, and 3.2 the baryon fraction was obtained 
as described in the final part of Sect. 2. Ordered angular 
momentum is fixed to typical values obtained by Gottlober 
& Yepes (2007), and random angular momentum is used as 
a parameter. 

The results of our calculation are plotted in Fig. la-c. 
The ordered angular momentum used to obtain the solid 
line in Fig. la-lc (obtained by the tidal torque theory as 
in DP09) is characterized by A = 0.03[f]. The baryonic frac- 
tion is F B , = M b /M 500 = 0.15, or f d = F b /f b ~ 0.88 (see 
McGaugh et al. 2010), and the typical random angular mo- 
mentum is j*. The solid line in Fig. la-c is obtained using 
these typical values. The solid line, obtained as described is 
chosen as the reference halo. 

Differently from DP12, in the present paper, we will 

5 Note that even if we studies haloes of mass in the range 
10 14 — 1O 15 M0, and this imply differences in ordered angular 
momentum, the spin parameter, A has small differences. 



study the change of the inner density profile slope just when 
the random angular momentum changes, since differently 
from dwarf and spiral galaxies, clusters of galaxies are non- 
rotationally supported (the mean rotational velocity, and 
ordered angular momentum of most luminous clusters are 
much less than their velocity dispersions). 

In fig. la, we fix the value of fd to the value previously 
discussed, and vary the amplitude of the random angular 
momentum, j. The short-dashed (dotted) line in Fig. la, rep- 
resents the case j = J./1.5 (j = j*/2), where is the value 
of j for the reference case (solid line). As expected, the den- 
sity profile steepens when j decreases (a ~ 0.6 (solid line), 
0.75 (short-dashed line), 0.95 (dotted line)). This steepen- 
ing of the profile with j agrees with several previous re- 
sults (Sikivie et al. 1997; Avila-Reese et al. 1998; Nusser 
2001; Hiotelis 2002; Le Delliou & Henriksen 2003; Ascasi- 
bar, Yepes & Gottleber 2004; Williams et al. 2004). The 
long-dashed (dot-dashed) line in Fig. la, represents the case 
j = j* x 1.5 (j = j» x 2). As expected, larger values of j 
produce flatter profiles (a ~ 0.4, 0.2, for the long-dashed, 
and dot-dashed line, respectively). 

Fig. lb, is the same as Fig. la, but the halo mass is in 
this case 10 15 Mq. Short-dashed (dotted) line is character- 
ized by a ~ 0.93 (1.05), while long-dashed (dot-dashed) line 
is characterized by a ~ 0.52 (0.32). 

A comparison of Fig. la and Fig. lb, shows that the 
1O 14 M0 halo is less steep than the 1O 15 M0 one. The steep- 
ening of the cusp with mass is in agreement with previous 
studies (Cen & Ostriker 2000; Ricotti 2003). Both Fig. la 
and Fig. lb, show that haloes with higher angular momen- 
tum have smaller inner slopes. The steepening of the profile 
with mass is due to the fact that higher density peaks (larger 
43) are usually progenitors of more massive haloefl which 
are characterized by a larger central density contrast. As a 
consequence, a generic shell inside this peak will fill more 
strongly the central potential and will expand less than if it 
was located in a smaller density peak. The angular momen- 
tum acquired during expansion will be reduced and haloes 
will be more concentrated. It is important to notice that the 
quoted trend of increased central concentration as a func- 
tion of mass applies only to halos that started out as peaks 
in the density field smoothed with a fixed Rf scale (see Ap- 
pendix B of DP09). Our conclusions do not mean that, for 
example, clusters of galaxies will be more centrally concen- 
trated than galaxies, since different smoothing scales would 
apply in the two cases. The angular momentum dependence 
of the slope arises from the fact that less massive objects 
are originated from peaks having a smaller height. They 

6 v = <5(0)/<t, where <5 is the mean fractional density excess inside 
a shell, and a is the mass variance filtered on a scale Rf (see 
Appendix B of DP09). 

7 As discussed by Peacock & Heavens 1990; Del Popolo & Gam- 
bera 1996 or Gao & White (2007) (Fig.l), modelling the peaks as 

triaxial spheroids one obtains a peak mass M = 2 „ ( 47r//3 ) p k S* 

for 0.5 ^ 7 ^ 0.8 (where 7 and ij* are spectral parameters given 
by Eq. (B6 and B7) in DP09). The peak mass is an increasing 
function of v. It is reasonable that lower v peaks should have 
lower mass; peaks with v ~ will tend to sit in regions of larger 
scale undcrdensity (cancelling the small-scale overdensity ) , and 
hence the ~ pcR\ of material, which initially surrounds the peak, 
may not be accreted following central collapse. 
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Figure 1. Changes of the density profile with baryonic fraction and random angular momentum j. Panel (a) refers to density profiles 
with mass 10 14 Mq varying only j. The solid line is characterized by Fg t = M^/M^oo = 0.15 (fd, = 0.88), and j = j t . The upper 
short-dashed line, and the dotted line are characterized by f d = f dt , and j = j',/1.5, and f d = / dt , j = j»/2, respectively. The long 
dashed line, and dot-dashed line, are characterized by fd = fd,, and j = j* X 1.5, and fd = fd,, j = j, X 2, respectively. Panel (b), 
represents the density profile for haloes of mass 10 15 Mq. Panel (c) plots density profiles for haloes of mass 10 14 Mq varying both j and 
baryonic fraction. The solid line, short-dashed line, dotted line, are characterized by fd = fd,, j = j,; fd = /d„/1.5, j = j»; fd = /d,/l-5, 
j = j*/2, respectively. The long-dashed line, dot-dashed line, are characterized by fd = 1.5 X fd t , j = j*; fd = 1-5 X fd t , j = j* X 2, 
respectively. 



acquire less ordered, h, and random, j, angular momenta 
(Ryden & Gunn 1987). The larger is the angular momen- 
tum, the more particles constituting the shell will remain 
closer to the maximum radius, producing a flatter density 
profile. Summarizing, particles having larger h and j tend 
to live far from the center, so not contributing to the central 
density. At the same time, a larger value of mass produces 
an increase of the central density contrast, and a steeper 
profile. 

The baryon content of the halo is another parameter 
playing a key role in shaping density profiles. In Fig. lc, 
we studied how changes in the baryonic fraction change the 



profile of a 10 Mq halo. Similar results, not plotted, are 
valid for 10 15 M Q haloes. 

In Fig. lc, the solid line is the same as Fig. la, while the 
short-dashed (dotted) line represents the case fd = /d»/T5 
and j = j* (f d = fdJl-5 and j = j„/2), where f d , = 
F B Jh = 0.15/0.17 = 0.88. The long-dashed (dot-dashed) 
line represents the case fd — fd, x 1-5 and j = j, (fd = 
fd, x 1.5 and j = j, x 2). 

The short-dashed line shows that a decrease of the bary- 
onic fraction, from f d , to / d „/1.5 produces a steepening of 
the density profile, to a ~ 0.7. If also j is decreased, to j*/2, 
a further steepening of the profile is produced (a ~ 1.1). In 
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Figure 2. Changes of the inner slope a in terms of the ratio of barionyc to total mass, M^/M^qq. Dots with error- bars refer to the case 
Id = fd, > j = j*- The triangles with 1 <r error- bars refers to the case f d = f dt , j = j,/2, and the squares with error- bars to the case 
fd = fd, > J = J* x 2. 



the reverse case, the long-dashed line, representing the case 
fd = fd, x 1.5 and j — j*, the profile flattens (a ~ 0.35), 
and a further flattening of the profile is produced (a ~ 0.1 if 
also, j is increased to j = j» x 2). The clear tendency to have 
flatter inner profiles for larger baryon content is connected 
to the fact that the presence of a larger quantity of baryons 
guarantees a larger transfer of energy and angular momen- 
tum from baryons to DM with the results that DM moves 
to larger orbits reducing the inner density. It is important to 
notice, that an increase in angular momentum has a more 
important effect on the inner density profile slope than an 
increase in the baryon fraction (an increase of a factor of 2 
in random angular momentum has a larger effect than an a 
similar increase in the baryon fraction). 

In the case that the system does not contain baryons, 
like dissipationless N-body simulations, the density profile 
converge to the Einasto profile (as shown in Fig. Id of 
DP09), characterized by a logarithmic slope that varies con- 
tinuously with radius (Navarro et al. 2004, 2010). 

In Fig. 2, we plot how the negative logarithmic slope 
of the DM density at small radii (10 _2 ri,i r ), c 



d log p 
d log r ' 



varies with the baryonic content in haloes, Mfj/Msoo- The 
dots with 1 a error-bars represent the relation a-Mb/Msoo 
in the case fd = fd, an d j = j*. The plot shows a de- 
crease of a with increasing value of the baryonic fraction, 
fd, as already noticed in Fig. 1: the final density profiles are 
shallower when F;, (or similarly f d ) is higher. The triangles 
(squares) with 1 a error-bars, represent the same relation in 
the case f d = f d , and j = j,/2 (f d = f d , and j = j, x 2). 
The plot shows that for a given f d one can obtain different 
slopes according to the value of random angular momentum 
of the system, and this suggests that the cluster mass and 



baryonic fraction are not the only fundamental parameters 
in shaping the density profile. A fundamental role is played 
by the orbital properties of the objects constituting the clus- 
ter. Galaxies moving on different orbits heat the DM at dif- 
ferent rates. Galaxies having larger kinetic energy transfer, 
through dynamical friction, a larger amount of energy to 
DM. So, if baryon content has a certain importance in de- 
termining the final DM distribution, the orbital parameters 
of galaxies and their dynamics has a similar or even a larger 
importance. In a recent paper, Laporte et al. (2012) studied 
how the evolution of galaxies influence clusters of galaxies 
formation. They followed galaxies evolution through colli- 
sionless mergers and found that cusps with 0.8 < a < 1.3 
can be flattened to 0.3 < ot < 0.9 at the innermost resolved 
radii. This result is in qualitative agreement with the results 
of our paper. However, in their study, the density profiles de- 
pend not only on the orbital properties of galaxies, as in the 
present study, but also on the internal structure of galax- 
ies. Less tightly bound galaxies (i.e., less compact) are more 
easily stripped with respect to more tightly bound ones, and 
as a result the more is the quantity of less compact galaxies 
in the cluster the more quantity of stellar mass is stripped 
and deposited on the BCG. The effect on the slope of the 
DM profiles is that clusters containing more compact galax- 
ies have shallower inner profiles and vice versa. Even if the 
model is an improvement on that of El-Zant et al. (2004), 
it has several limitations as recognized by the authors, and 
certain aspects of matter distribution in their simulations 
are unrealistic. This is one of the reasons they did not even 
attempt to reproduce observations of clusters. Moreover, the 
flattening of the inner slope of the density profile of clusters 
is not only due to the effect of galaxy evolution. Dissipa- 
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tional processes have an important role in shaping galaxies 
at early times. Baryons are more important than DM in the 
inner parts of galaxy clusters and alter the DM profile in 
many ways, as described in the previous part of the present 
paper. 

In Fig. 1, and Fig. 2, we discussed how angular mo- 
mentum shapes the density profiles, and that haloes charac- 
terized by a larger angular momentum have flatter density 
profiles. This result could erroneously lead us to think that 
angular momentum is the primary effect driving the shap- 
ing of the density profiles, and the principal reason giving 
rise to the difference between density profiles found in sim- 
ulations and those found in the present paper. This is not 
the case. We have to recall that even if we are studying 
the effect of angular momentum change on the inner den- 
sity slope, our system is composed of baryons and DM, and 
as a consequence our results are different from dissipation- 
less N-body simulations. In our model, angular momentum 
is not working alone in shaping the density profiles, dynam- 
ical friction transfer angular momentum from baryons to 
DM. The larger is the angular momentum in the system, 
the larger is the quantity that can be transferred to DM 
with a consequent larger flattening of the density profile. 
I want also to add that the flattening of density profiles 
with increasing magnitude of angular momentum was also 
obtained in several studies even of systems not containing 
baryons (Avila-Reese et al. 1998, 2001; Subramanian et al. 
2000; Nusser 2001; Hiotelis 2002; Le Delliou & Henriksen 
2003; Ascasibar, Yepes & Gottleber 2004; Williams et al. 
2004; Ascasibar, Hoffman & Gottlober 2007). In peculiar, 
Ascasibar, Yepes & Gottleber (2004), compared the analyt- 
ical profiles of SIM which included non-radial motions with 
a set of high-resolution N-body simulations, showing that 
angular momentum is responsible for the shape of the den- 
sity profile near the centre. Ascasibar, Hoffman & Gottlober 
(2007) used N-body simulations to show that SIM gives den- 
sity profile in good agreement with N-body simulations and 
that larger values of angular momentum produce a flatten- 
ing of the profile. 

Angular momentum acts also on galaxies, and if it was 
the only effect taken into account in our model, it would 
delay their collapse, keeping them away from the center. In 
reality galaxies, as baryons, are also subject to dynamical 
friction which makes them loose their energy and sink to 
the center (Kashlinsky 1987). 

The reason why our model gives different results from 
dissipationless N-body simulations was already discussed in 
DP09, and it is not connected to the way our model takes ac- 
count of angular momentum. The issue is not if our model 
takes into account angular momentum more reliably than 
N-body simulations, but the difference among the system 
described in our model and that of dissipationless simula- 
tions. In DP09 (beginning of page 2101), we described the 
two-fold effect of baryons presence, namely the steepening 
of the density profile due to AC and the flattening of the 
inner density profile due to the transfer of angular momen- 
tum from the baryons to the DM. In collisionless N-body 
simulations, in which baryons are not present, this compli- 
cated interplay between different effects is not taken into 
account. Only recent SPH simulations (Romano-Diaz et al. 
2008; Governato et al. 2010) took account of the baryonic 
component, finding the quoted cusp erasing. 



Finally, we recall that a check of the model against the 
Governato et al. (2010) SPH N-body simulations, has been 
performed in DP12 (Fig. 3). 

3.2 Comparison with observed clusters of galaxies 

At this point, we compare the results of our model with 
density and mass profiles of observed clusters, ranging from a 
fiat inner profile, namely that of A611, to intermediate ones, 
namely A383, to steep profiles, namely MACSJ1423. 8+2404 
and RXJ1347.5-1145. 

A611, A383, and RX J1133 are three of the clusters 
studied in a series of papers of Sand (S02, S04, and S08) and 
Newman (N09, N011). The quoted papers, followed a pro- 
posal of Miralda-Escude (1995), evidencing that if for a clus- 
ter we can add to the information on tangential and radial 
arc the information concerning the potential of the central 
galaxy (BCG), stronger constraints on the inner profiles of 
DM and luminous components can be obtained. S02 studied 
MS 2137-23, while S04 six clusters (MS2137-23, A383, A963, 
RX J1133, MACS 1206, and A1201), and S08 improved the 
results concerning MS2137-23 and A383. Nil improved the 
constraints previously obtained in S08 concerning A383, and 
N09 studied A611. 



3.2.1 A611 

A611 has been studied in several papers (OWLS team; Bona- 
mente et al. 2004; Schmidt & Allen 2007; Hurley- Walker 
et al. 2011 (AMI consortium)). N09, combined weak lens- 
ing from multicolor Subaru imaging, strong lensing (Hub- 
ble Space Telescope), and stellar velocity dispersion mea- 
sures (Keck Telescope), sampling the dark matter profile 
from ~ 3 kpc to 3.25 Mpc. In order to accurately com- 
pare observational data and N-body simulations, one needs 
a) to separate DM from baryons, which, nevertheless they 
are a small part of the total mass of clusters, are usually 
dominating in the inner kpc scale; b) to probe mass on all 
scales. To this aim it is necessary to consider a wide dynamic 
range in radius. Concerning the mass distribution in the 
range 150 kpc- 3.25 Mpc, N09 found values of r s = 320±i 4 g 
kpc, and c = 5.1 jl x \ (in agreement with Schmidt & Allen 
2007). The distribution of the stellar mass of BCG, mea- 
sured through strong lensing and stellar velocity disper- 
sions, is well fitted by a ii 1 / 4 law in the 3-20 inner kpc. 
The BCG is characterized by L B = (5.6 ± 0.8) x 10 11 M Q , 
and M+.bcg/Lb = 2.7t° g (N09). Combining weak lens- 
ing, strong lensing, and stellar kinematics data, then they 
concluded that a < 0.3, (< 0.56, < 0.65) at 68%, (95%, 
99%) CL. Schmidt & Allen (2007) found a value of the slope 
a = 0.64+ ' 94 . An important result obtained by N09 is that 
using a subset of their data coming from just one technique 
(e.g., weak lensing, or strong lensing), the data are well fit- 
ted by a NFW model, but NFW model is unable to fit all 
data simultaneously. This imply the need to use several com- 
bined techniques and large radial scales. We compared the 
result of our model to N09 DM density profile in Fig. 3a. The 
cluster mass is estimated as M„, r = 6.18ti'gi x lO 14 /i _1 M0 
(Schmidt & Allen 2007), M 200 = S-^ttl x 1O 14 /i _1 M (Ro- 
mano et al. 2010). We recall, that, as previously reported, 
all masses in the paper were converted to Msoo- In order to 
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Figure 3. Fig. 3a plots the density profile of A611. The grey solid line is the NFW profile, the red lines bracket the 68% confidence 
region of the density profile obtained by N09, while the black band is the result of our model. Fig. 3b represents the mass profile of A383. 
The blue band is the observational result of Nil, while the black band the mass distribution obtained with our model. Fig. 3c represents 
the mass profile of MACS J1423. The azure band with errorbars is the observational result of Morandi, Pedersen & Limousin (2010), 
while the black band is the result obtained with our model. Fig. 3d represents the density profile of RX J1133. The solid curve and the 
shaded region are the S04 observational result. The black dashed band is the result of our model. The confidence regions in Fig. 3a-c are 
1 a while those in Fig. 3d are 2 a. 



determine the DM density profile, we used our model. The 
baryonic fraction obtained as described in the final part of 
Sect. 2, gives a value of M b /M 500 ~ 0.15 (f d ~ 0.88). The 
value of the baryonic fraction is larger with respect to that 
of RX J1133, and MACS J1423 (as we shall see in the fol- 
lowing), and this is one reason why one should expect a 
flatter inner profile with respect to the two quoted clusters, 
well fitted by a NFW model. In Fig. 3a, the grey line is the 
NFW fit, the red lines bracket the 68% confidence region of 
the density profile obtained by N09, while the black band is 
the result of our model (68% CL). Looking at Fig. lb, for a 
cluster having a mass ~ 10 15 M Q and fd = fd, , and j = j», 
the inner slope is of the order of 0.7. Taking account the 
fact that A611 has a mass smaller than 1O 15 M0, one would 
obtain a value of a ~ 0.65. In order to obtain the good fit 
plotted in Fig. 3a, we had to increase the magnitude of j to 
2 x j*. As shown in Fig. 1, and also in Williams et al. (2004) 
(Fig. 6), increasing the value of j produces a flattening of 
the profile. 



3.2.2 A383 

A383, is one of the most X-ray luminous clusters in the red- 
shift slice 0.17 < 2 < 0.26, at z = 0.189. It is characterized 
by numerous strong lensing features and a dominant cen- 
tral galaxy. An analysis of the mass distribution in the core, 
built up through a lens model of the cluster, was carried 
out by Smith et al. (2001, 2005). The inner density profile, 
was evaluated by them to be a = 1.3 ± 0.04. The cluster 
was one of the quoted six clusters studied in S04, combin- 
ing results from stellar velocity dispersion data, and lensing. 
The cluster has both a radial and tangential arc. The dis- 
persion velocity profile, and the density profile are plotted 
in Figs. 7, 8 of S04. The inner DM slope obtained by S04 
was a = 0.38±g:g|(±g:i|) (68% (95%) CL). The results of 
S04, and consequently the constraints on A383 inner slope, 
were criticized by Meneghetti et al. (20070, on the base 



8 Notice that Meneghetti uses traxial clusters similar to those 
obtained in simulations, while those in S04 are more spherical. 
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that S04 used axially symmetric models in order to obtain a 
description of the six clusters mass model, while a pseudo- 
elliptical lens model (including ellipticity in the lensing po- 
tential) gives rise to steeper density profiles. The previous 
criticisms were tackled in S08, in the case of A383, finding 
again a flat inner DM slope (a = 0.45lg 25)- 

A further improvement to the model was presented in 
Nil, in which the authors, used the same technique used 
in N09. They also combined X-ray and lensing constraints, 
in order to measure the DM elongation along the line of 
sight. In this way, it was further possible to tackle an im- 
portant systematic uncertainty in the mass profile determi- 
nation (see Morandi et al. 2010). As a result, they further 
improved the S08 constraints, who has not taken into ac- 
count triaxiality. They obtained a value of a < 0.70 (68% 
confidence), a < 1 (95% confidence), and a best fit (inferred 
from weak, strong lensing, kinematics, and X-ray data) of 
a = 0.5910 35 J- n order to compare the mass profile ob- 
tained by Nil, with our model, we need mass and baryonic 
content of the cluster. Zitrin et al. (2011) and Schmidt et al. 
(2007) found values of M vir = 6.26±o;25 x 10 14 /i _1 M Q , and 
M vir = 6.87±i;gg x 10 14 /i -1 M s , respectively. The mass of 
the BCG, within ~ 19 kpc, is 1.14 ± 0.03 x 10 12 Af o (Zitrin 
et al. 2011) in agreement with Nil. The baryonic fraction 
was obtained as for A611. The value of M&/M500 is — 0.15 
(fd ~ 0.88). In Fig. 3b, we compare the mass distribution 
obtained with our model (black band), with that of Nil 
(blue band). The good fit to the mass distribution of A383 
is obtained with the typical value of random angular mo- 
mentum. This shows that even if A383 has a mass and bary- 
onic fraction close to that of A611, the smaller value of j give 
rise to the steeper profile observed (a ~ 0.6). 

3.2.3 MA CSJ1 423.8+2404 

The cluster of galaxies MACS J1423. 8+2404 (in the fol- 
lowing referred as MACS J1423), is a triaxial LJ massive 
(M = 4.52io.64 x 1O 14 M , Schmidt & Allen 2007) strong 
cooling core sourc (0 (Morandi et. al. 2007b) at z = 0.539. It 
is very relaxed (Kartaltepe 2008) with very low central tem- 
perature (~ 2 keV). Morandi, Pedersen & Limousin (2010) 
performed a joint analysis (X-ray, strong, and weak lensing) 
to determine DM density profile and ICM parameters. The 
BCG mass is 5 x lO n M0. The three-dimensional analysis 
of the cluster, lead (Morandi, Pedersen & Limousin 2010) 
to measure an inner slope a = 0.94 ± 0.09 smaller than 
the two-dimensional spherical modeling, giving a value of 
a = 1.24 ± 0.07. Similarly to the previous clusters, we used 
our model to determine the mass profile of the cluster. The 
value of the baryonic fraction, inferred as in A383, and A611, 
gives Fb — 0.1 (fd — 0.59), and the mass profile, plotted in 
Fig. 3d (black band), in good agreement with Morandi, Ped- 
ersen & Limousin (2010) (azure band with errorbars), was 
obtained assuming j — j*/1.5. The effect of steepening of 

9 The authors are also studying other nine clusters with the same 
level of precision. 

10 DM halo axial ratios are 1.53 ± 0.15, on the plane of the sky, 
and 1.44 ± 0.07, along the line of sight (Morandi. Pedersen & 
Limousin 2010). 

11 i.e., the central cooling time is much smaller than the age of 
universe. 



the profile with decreasing j was pointed out by Williams et 
al. (2004) (their Fig. 6). Reducing the angular momentum 
of dark matter particles, reduces random velocities which 
results in steeper central density slopes. 

3.2.4 RX J113S 

RX J1133 accurately studied in S04, is a cluster with total 
mass ~ 3.2 x 1O 14 M0 (Cardone, Piedipalummbo, & Tortora 
2005; Lackner & Ostriker 2010 (private communication)), 
and a BCG with mass 3 x 10 11 Mq (Lackner & Ostriker 
2010 (private communication)). S04 identified a tangential 
and radial arc and measured the velocity dispersion profile 
of the BCG. S04 disentangled the DM and baryonic com- 
ponent, finding density profiles of DM, luminous matter, 
and total matter (see S04, Fig. 8). The inner slope of the 
DM profile was found consistent with that of a NFW profile 
(a = 0.99^0.14 )■ In Fig. 3d, we plot the DM density profile 
of RX J1133. Fig. 3d represents the density profile of RX 
J1133. The solid curve and the shaded region is the S04 ob- 
servational result. The black dashed band is the result of our 
model, using the previously given mass for the cluster, and 
the baryonic fraction calculated as for the previous three 
clusters. The obtained value Fb — 0.1 is smaller than that of 
A383 and A611, implying a steeper slope with respect to the 
two quoted clusters. The model of this paper gives a very 
good fit to the DM density profile of RX J1133, obtained 
by S04, if j — j*/2. The reduction of random angular mo- 
mentum is justified as in the case of MACS J1423. Finally, 
we notice how the result shows that a NFW model is also a 
good fit to RX J1133 density profile. 

The confidence regions in Fig. 3a-c are 1 a while those 
in Fig. 3d are 2 a. 



4 DISCUSSION 

4.1 What are the origins of the differences? 

From the results of the previous sections, we arrive to two 
important conclusions, that we will discuss: a) not all clus- 
ters density profiles are fitted by the NFW model. Some of 
them have flat inner density profiles (e.g., A611), some have 
intermediate slopes, between pseudo-isothermal profiles and 
NFW profiles, (e.g., A383), and others are well fitted by 
NFW profiles (e.g., RX J1133); b) the ACDM model is un- 
able to describe the density profile of some clusters. 

In connection with the first conclusion, a connected 
question is the following: what is the cause of the differ- 
ences of density profiles in clusters? In the study of the role 
of baryons and random angular momentum developed in 
Sect. 3.1, we saw that the causes that give rise to differ- 
ent slopes in clusters are three, namely the mass (clusters 
with larger mass have steeper profiles), the baryon content 
(clusters richer in baryons have flatter profiles), and random 
angular momentum, j (larger values of j implies flatter pro- 
files). We also saw that random angular momentum have a 
stronger role than baryons in shaping density profiles. A383 
and A611 have a similar mass (M vir ~ 6 x lO 14 /t _1 M ), 
and a similar baryonic fraction (Fb — 0.15 (fd — 0.88)). 

12 Obviously issue a. and b are strictly connected. 
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RX J11333 and MACS J1423 have slightly smaller masses 
(~ 3 — 4 x 10 14 h~ 1 M@), and smaller baryonic fraction 
(Fb — 0.1 (fd — 0.59)). The difference in baryonic frac- 
tion is one of the cause of the steeper profile of RX J11333 
and MACS J1423 with respect with the other two clusters, 
in agreement with observations. We should also not forget 
that, similarly to dwarf galaxies, in order two clusters have 
a similar density profile, it is not only important that the 
baryonic fraction is similar for the two objects, but the way 
baryons are distributed has also an important role. 

In particular, the larger or smaller concentration of gas 
and stars in the central ~ 10 kpc, in the BCG, is also impor- 
tant. As shown by Schmidt & Allen (2007) (Fig. 4), the inner 
slope of clusters of galaxies (MS 2137.3-2353, in particular) 
decreases with increasing values of the central mass. In the 
inner ~ 10 kpc, A611 has a Mbcg /M v i r ~ 0.004, using the 
data of N09, and of Schmidt & Allen (2007), concerning 
Mbcg /M vir . A383 has M B cG/M vir ~ 0.002 (within ~ 19 
kpc), using Nil; Schmidt & Allen (2007), and Zitrin et al. 
(2011), data. RX J1133 has M B CG/M vir ~ 0.001 (Cardone, 
Piedipalummbo, & Tortora 2005; Lackner & Ostriker 2010 
(private communication)), and MACS J 1423 has a similar 
value. The larger value of ratio of the inner stellar mass to 
total mass of A611 with respect to A383, and the larger 
ratio of A611, and A 383 with respect to RX J1133 also 
implies that the inner slope of A611 is flatter than that 
of A383, and flatter than those of MACS J1423, and RX 
J1133. 

It is interesting to discuss more in detail the role of the 
total baryonic mass, Mf,, and the central one. As shown in 
El-Zant et al. (2001, 2004), Romano-Diaz et al. 2008, DP09, 
Governato et al. (2010), baryons presence produces in gen- 
eral a flattening of the density profile. The baryon diffused 
component produces rounder halos and less triaxial halos 
than those seen in DM simulations (Gustafsonn et al. 2006; 
Debattista et al. 2008; Abadi et al. 2010), and this change of 
shape has also influence on dynamical mass estimate (N09). 
The final configuration of a cluster is fixed by the initial 
quantity of baryons present in the proto cluster and by col- 
lapse/formation process. 

From this considerations, we should expect, in hierar- 
chical formation models, that the final central baryonic con- 
tent and the BCG mass is somehow correlated with baryonic 
and total cluster mass. This is in fact the shown by 

Whiley et al. (2008), who using the models of de Lucia & 
Blaizot (2007), found that Mbcg oc M°( 4 or M°f J depend- 
ing on the feedback model used. Also a correlation between 
BCG luminosity and cluster X-ray luminosity was found by 
several authors (Schombert 1988; Edge 1991; Edge & Stew- 
art 1991; Hudson & Ebeling 1997). Whiley et al. (2008) 
measured the quoted correlation as Mbcg oc j^f0- 12 ± - 03 f or 
K band magnitudes inside a diameter of 37 kpc (radius of 
13/i _1 kpc). Brough et al. (2008) found Lbcg oc M°' 11±0 ' 10 
at K band inside 12/i _1 kpc (several other results are given 
in Lin & Mohr 2004; Popesso et al. 2007; Yang et al. 2008; 
Haarsma et al. 2010). 

In Fig. 4a, we calculated the baryonic mass in the inner 
10 kpc and plotted it against the total mass. The dots with 
1 a errorbars represent the quoted relation for fd = fd, and 
j = j,. The triangles with errorbars represent the quoted 
relation for fd = fd, and j = j, x 2. The squares with 



errorbars represent the quoted relation for fd = fd, and 
j = j,/2. The plot shows a similar correlation as that found 
in Whiley et al. (2008) 

Mb ' m ~ i 3x in 6 ( M50 ° V 4 (17) 

10" M - LJX 10 {w^MqJ {17) 

Fig. 4b, plots the inner slope (defined as in Fig. 2) in 
terms of the ratio Mj, iin /Af 5 oo. Dots with 1 a error-bars rep- 
resent the quoted relation in the case fd = fd, and j = j, . 

In Fig. 4b, it is also important to note that the funda- 
mental parameter influencing the inner slope of the density 
profile is Mi,^ n /M v i r and not just the baryonic mass Mj, jin . 
In fact, if two clusters have the same values of Mb,i n , but 
different M V i r , the one having smaller value of M„; r would 
have a flatter profile, since the role of the central baryonic 
component is larger. This was also noticed by Schmidt & 
Allen (2007) (Fig. 4), in the case of MS2137.3-2353, which 
is one of the least massive clusters in the sample that they 
used. They showed that increasing the values of the central 
stellar mass produces larger effects on the inner slope value, 
with respect to other cluster, because, as reported, the to- 
tal mass of MS2137. 3-2353 is smaller than the other sample 
components. 

The other quantity of fundamental importance in shap- 
ing the density profiles is the random angular momentum. 
We have seen that a good fit to the mass/density profile 
of the studied clusters implies different values of j. Fig. 2 
clearly shows that cluster mass and baryonic fraction are 
not the only fundamental parameters in building up the clus- 
ter structure. Random angular momentum, j, strongly influ- 
ences cluster formation. During the collapse, particle follows 
orbits connected to the value of j that they have, and par- 
ticles endowed with larger kinetic energy will transfer more 
energy (through dynamical friction) to DM with the results 
that DM will expand reducing the inner density. So, as pre- 
viously reported, galaxies orbital parameters, and dynamics 
have a fundamental role in density profile formation, even 
larger than the baryonic content. 

4.2 Problems for the ACDM 

In the previous subsection, as in the introduction, we pointed 
out that ACDM model is unable to explain the flat density 
profiles observed in dwarf galaxies, and to describe the den- 
sity profile of some clusters, having too flat inner slopes with 
respect to the ACDM model predictions. However this prob- 
lem does not necessarily imply a problem for CDM model 
(e.g. DP09; Governato et al. 2010). In the clusters that we 
studied, only RX J1133 is in agreement with ACDM model 
predictions, while the inner slope of the density profiles of 
A611 and A383 contradicts the N-body simulations results. 
Although more recent simulations (e.g. Navarro et al. 2010) 
found density profiles characterized by a continuous flatten- 
ing of the inner slope, this does not solve the contradic- 
tion between observations and simulations, since the quoted 
clusters have smaller slopes than the minimum value of the 
slope (a = 0.8 at 120 pc) found in high-resolution dissipa- 
tionless N-body simulations (e.g., Stadel et al. 2009). The 
quoted discrepancy is not restricted to the clusters stud- 
ied in the present paper. On the observational side, as dis- 
cussed in introduction, X-ray observations (Ettori et al. 
2002; Arabadjis, Bautz & Garmire 2002; Lewis, Buote & 
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Figure 4. Fig. 4a: baryonic mass in the inner 10 kpc in terms of the total mass. Dots with 1 a error-bars, triangles with errorbars, and 
squares with error-bars represent the quoted relation in the cases f d = and j = j„; f d = f dt and j = X 2; and f d = f dt , and 
j = j'»/2, respectively. Fig. 4b: inner slope in terms of the ratio M^i n j M^qq. Dots with error-bars represent the quoted relation in the 
case f d = f dt and j = j, . 



Stocke 2003), lensing (Tyson et al. 1998; Smith et al. 2001; 
Dahle, Hannestad & Sommer-Larsen 2003; S02; Gavazzi et 
al. 2003; Gavazzi 2005; S04; Bradac et al. 2008; Limousin et 
al. 2008), dynamics (Kelson et al. 2002; Biviano & Salucci 
2006), or studies combining several techniques (e.g. S02; S04; 
S08; N09, Nil), lead to large scatter in the value of a from 
one cluster to another. If some scatter can be explained, as 
reported in introduction by different/limited dynamic range 
in radius in different studies, BCG role not taken into ac- 
count, or as pointed out by Morandi, Pedersen & Limousin 
(2010), to standard (simplified) spherical modeling of clus- 
ters, as well as the degeneracy of a with parameters like c 
and r s , it is difficult to explain the large slope difference 
between some clusters (e.g., A611). If the scatter in the in- 
ner slope of density profile is not merely caused by limits in 
the techniques used, some other reason causes it. We could 
conclude that either some assumptions of ACDM model are 
incorrect, or given the several and noteworthy evidences sup- 
porting ACDM on large scales, another possibility, moreover 
studied in DP09, is that baryonic physics is a fundamental 
issue in clusters formation. Probably, as noticed in DP09, 
the quoted discrepancy is merely due to the fact that we are 
comparing two totally different systems: the one generated 
by dissipationless simulations, not including baryons, with 
real structures whose physics is not just the dissipationless 
physics typical of DM. 

To start with, we would recall that dissipationless simu- 
lations does not include baryons, which are usually dominant 
in the inner part of clusters 13 !, and its presence strongly in- 
fluence the DM distribution. If stars form earlier than DM, 
baryons will compress the DM (adiabatic contraction), giv- 
ing rise to steeper profiles (Blumenthal. et al. 1986; Gnedin 
et al. 2004; Gustafsonn et al. 2006). However, adiabatic con- 
traction can be counteracted by heating of DM due to dy- 
namical friction with cluster galaxies (El-Zant et al. 2001, 

13 For example, within 10 kpc, the total density distribution is 
dominated by the BCG in A383, RX J1133, A1201, A963, MACS 
1206, MS2137-23 (S04). 



2004; Nipoti et al. 2003; Romano-Diaz et al. 2008). Several 
studies have tried to rescue the ACDM paradigm without 
drastic changes to the ACDM physics, and studying baryon 
physics and stellar processes in the inner parts of galax- 
ies and clusters. Interactions of the DM with a stellar bar 
(Weinberg & Katz 2002; McMillan & Dehnen 2005), decay 
of binary black hole orbits after galaxies merge (Milosavl- 
jevic & Merritt 2001), baryon energy feedback from active 
galactic nucleus (Peirani et al. 2008), the already quoted 
dynamical friction of stellar/DM clumps against the back- 
ground DM halo (El-Zant et al. 2001, 2004; Romano-Diaz 
et al. 2008, 2009), random bulk motions of gas in primordial 
galaxies, driven by supernova explosions (Mashchenko et al. 
2006), removal of low-angular momentum gas (Governato et 
al. 2010), are some of the solutions proposed. Moreover, Zap- 
pacosta et al. (2006) concluded, through X-ray observations 
of A2589, that processes in clusters of galaxies counteract 
adiabatic contraction. Infalling DM subhaloes, according to 
their mass, could steepen or flatten the DM cusp (Ma & 
Boylan-Kolchin 2004), and Nipoti et al. (2004) got similar 
results according to the baryon fraction of infalling galaxies. 

Another important caveat is that density profiles are 
not universal, and that in nature may exist a distribution 
of inner slopes broader than that predicted by dissipation- 
less N-body simulations. Every cluster has its own forma- 
tion/merger history leading to different cluster characteris- 
tics (see Navarro et al. 2010, and S04). Several studies argue 
against universality of DM profiles (Jing & Suto 2000; Sub- 
ramanian 2000; Ricotti 2003, 2004, 2007; Cen et al. 2004; 
Simon et al. 2005; Merrit et al. 2005; Graham et al. 2006; 
Schmidt et al. 2008; Del Popolo 2009; Ma et al. 2009; Host & 
Hansen 2010; Del Popolo 2012). N-body simulations as Jing 
& Suto (2000), Fukushige et al. (2003), Ricotti (2003), Ri- 
cotti & Wilkinson (2004) found inner slope variations from 
run-to-run, or with mass. Even Navarro et al. (2004, 2010) 
found dependence of the inner slope with mass, interpreted 
as a reflection of the trend between the concentration of 
a halo and its mass. Gao et al. (2008) in two very large 
cosmological simulations, found that density profiles devi- 
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ate slightly but systematically from the NFW form and are 
better approximated by an Einasto profile. Moreover, the 
shape parameter of the quoted profile changes with mass 
and redshift. Merrit et al. (2005, 2006) interpret the varia- 
tion in profile shape with halo mass, as an indication of the 
fact that ACDM halos have not a really universal profile, 
as already claimed by all the authors previously quoted. By 
using DP09, we studied, in Del Popolo (20f 1), the pseudo 
phase-space density, arguing against universality of density 
profiles constituted by dark matter and baryons. From the 
observational point of view, Simon et al. (2005), S02, S04, 
S07 studies, argue against universality of density profiles. 

If density profiles are not universal, is of fundamental 
importance to collect further high quality data, like those 
in A6II or A383, for larger samples of clusters (see N09, 
and Nil) to calculate not only the mean of the distribution 
but also their moments. At the same time, SPH simulations, 
similar to that of Governato et al. (2010) for dwarfs, could 
be run to study clusters formation and evolution, then com- 
paring the results with the measured distribution of inner 
slopes. 



5 CONCLUSIONS 

In the present paper we have studied how changes in ran- 
dom angular momentum and baryonic fraction affect density 
profiles of clusters with masses in the range 1O 14 _1O 15 M . 
The paper extends the study of DPI2, on dwarf galaxies, 
to cluster of galaxies. A reference density profile was calcu- 
lated using DP09, and then we studied how density profile 
changes when changing baryonic fraction and angular mo- 
mentum. Similarly to the case of dwarfs, the inner density 
profile steepens with increasing value of the halo mass, and 
with decreasing values of angular momentum and baryonic 
fraction (see Figs. I and 2). This is due to the fact that when 
more baryons are present the energy and angular momen- 
tum transfer from baryons to DM is larger, and DM moves 
on larger orbits reducing the inner density. Haloes consti- 
tuted only of DM have Einasto's density profiles. The pre- 
vious calculation was applied to four clusters (A6I f , A383, 
MACS JI423, and RX JII33). A6II has a flat inner pro- 
file, as found by N09, and a baryonic fraction Ft, ~ 0.15 
(fd — 0.88). Its density profile is re-obtained with a random 
angular momentum 2 x j* . A383 has a flat inner profile, but 
steeper than that of A6I I . It has a similar baryonic fraction 
and we re-obtain its mass distribution with the typical value 
j,. MACS J1423 and RX JI133 have steeper profiles, with 
RX JI133 well fitted by a NFW profile. They are charac- 
terized by a smaller baryonic fraction of A611 and A383, 
Ft — 0.1 (fd — 0.59) and their profiles are re-obtained with 
j — j*/1.5, and j = j*/2, respectively. Since the baryonic 
content in the inner kpcs of clusters can influence their inner 
slope, as described by Schmidt & Allen (2007), we studied 
how the baryonic mass in the inner ~ 10 kpc relates with 
the total mass of the cluster. We found a similar correlation 
to that found by Whiley et al. (2008), namely an increase of 
central baryonic mass with total cluster mass. Clusters hav- 
ing larger central baryonic mass have flatter profiles. We also 
found that a fundamental role is played by the orbital prop- 
erties of the objects constituting the cluster. So, if baryon 
content has a certain importance in determining the final 



DM distribution, the orbital parameters of galaxies, con- 
stituting the cluster, and their dynamics have a similar or 
even a larger importance. In summary, the density profile 
of clusters is strongly influenced by baryons, random angu- 
lar momentum, and from the orbital parameters of galaxies, 
constituting the cluster, and their dynamics. 

Differently from DPI2, in this paper we did not study 
the eventual correlation of density profiles with environment. 
In that paper, we found that environmental effects influence 
density profiles of dwarf galaxies. It would be interesting to 
deal with this issue in a future paper, even if this study is 
more complicated than in the case of dwarf galaxies. This 
study has a fundamental complication: in order to study an 
eventual environment- density profile correlation, we need a 
large sample of clusters from which to extract isolated and 
unisolated clusters. From an observational point of view, 
there have not been many studies on the effects of the en- 
vironment on the slopes of the dark matter haloes (but see 
DP12). Moreover, the concept of isolated or un-isolated is 
more or less well defined for dwarfs (e.g., Karachentsev et al. 
2004, Karachentsev & Kashibadze 2006), and some samples, 
even if not large, are already present in literature (WHISP; 
Swaters et al. 2002) and in the next few years some others 
will be added (the LITTLE THINGS and VL A- ANGST), 
for the case of clusters the situation is more complicated. 
To our knowledge, in literature Plionis et al. (2009) identi- 
fied sub-clusters within 2500 km/s of the main cluster, de- 
liberately excluding those cluster with distorted X-ray mor- 
phology and other indications of strong interactions. Hence 
their sample is conservative and we would expect them to be 
non-interacting in a significant way with neighboring ones 
(even if we cannot rule out that they are weakly interacting 
with a structure that isn't delineated due to sub-sampling of 
redshifts). From that, Pimbblet (2010) constructed another 
restricted, and more isolated sample. 

Clearly, much theoretical and observational work needs 
to be done in order to obtain insights into the role of the 
environment on the density profiles of cluster of galaxies. 



ACKNOWLEDGMENTS 

We would like to thank Tommaso Treu, Robert Schmidt, 
David Sand, Alister W. Graham, Claire Lackner, Andrea 
Morandi, Crescenzo Tortora, Igor Karachentsev, Valentina 
Karachentseva, Stefano Ettori, for stimulating discussions 
on the topics related to the subject of this paper. Finally we 
thank the referee, A. B. Newman for providing constructive 
comments and help in improving the contents of this paper. 



REFERENCES 

Aarseth S.J., Binney J., 1978, MNRAS 185, 227 

Abadi, M. G., Navarro, J. F., Fardal, M., Babul, A., & Steinmetz, 

M. 20010, MNRAS 407, 435-446 
Arabadjis, J. S., Bautz, M. W., Arabadjis, G. 2004, ApJ 617, 303 
Arabadjis, J. S., Bautz, M. W., Garmire, G. P. 2002, ApJ, 572, 

66 

Ascasibar Y., Yepcs G., Gottlobcr S., 2004, MNRAS, 352, 1109 
Ascasibar Y., Hoffman Y., Gottlober S., 2007, MNRAS, 376, 393 
Avila Reese V., Firmani C, Hernandez X., 1998, ApJ, 505, 37 



On the Density profile slope of Clusters of Galaxies 15 



Avila-Reese, V., Firmani, C, Klypin, A., & Kravtsov, A. 1999, 
MNRAS, 310, 527 

Avila-Reese, V., Colin, P., Valenzuela, O., DOnghia, E., & Fir- 
mani, C, 2001, ApJ, 559, 516 

Bardeen J.M., Bond J.R., Kaiser N., Szalay A.S., 1986, ApJ 304, 
15 

Barnes, J., & Efstathiou, G. 1987, ApJ, 319, 575 

Barrow, J.D., Silk, J., 1981, ApJ 250, 432 

Bartelmann M., Meneghetti M., 2004, A&A, 418, 413 

Bett, P., Eke, V., Frenk, C. S., Jenkins, A., Helly, J., Navarro, J., 

2007, MNRAS 376, 215 
Binney J., Silk J., 1979, MNRAS 188, 273 

Blais-Oucllettc S., Amram P., Carignan C, Swaters R., 2004, 

A&A, 420, 147 
Bonamente, et al. 2004, ApJ 614, 56-63 
Bradac, M., et al. 2005, A&A, 437, 49 
Bradac, M., et al. 2008, ApJ 681, 187 

Brough, S., Couch, W. J., Collins, C. A., Jarrett, T., Burke, D. J., 

& Mann, R. G. 2008, MNRAS, 385, L103 
Cardone, V. F., Piedipalumbo, E., Tortora, C, 2005A&A 429, 49 
Cardone, V. F., & Sereno, M. 2005, A&A, 438, 545 2003, ApJ, 

593, 26 

Ce n, R. Y., Dong, F., B ode, P., & Ostriker, J. P. 2004, 

|arXiv: astro-ph /0403352 1 
Dahle H., Hanncstad S., Sommcr-Larsen J., 2003, ApJ, 588, L73 
Debattista, V. P., Moore, B., Quinn, T., Kazantzidis, S., Maas, 

R., Mayer, L., Read, J., & Stadel, J. 2008, ApJ, 681, 1076 
de Blok W. J. G., Walter F., Brinks E., Trachternach C, Oh S-H., 

Kennicutt R. C, 2008, AJ, 136, 2648 
de Blok, W. J. C, & Bosnia, A. 2002, A&A, 385, 816 
de Blok, W. J. G., Bosnia, A., & McGaugh, S. 2003, MNRAS, 

340, 657 

Do Lucia, G. & Blaizot, J. 2007, MNRAS, 375, 2 
de Naray R. K., McGaugh S. S., de Blok W. J. G, 2008, ApJ, 
676, 920 

de Naray R. K., McGaugh S. S., Mihos J. C, 2009, ApJ, 692, 
1321 

Del Popolo, A., & Gambera, M. 1996, A&A, 308, 373 

Del Popolo, A., 2009, ApJ 698, 2093 

Del Popolo, A., 2010, MNRAS 408, 1808 

Del Popolo, A., 2011, JCAP 07, 014 

Del Popolo, A., 2012, MNRAS 419, 971 

Dressier A., 1978, ApJ 243, 2 

Edge, A. C. & Stewart, G. C. 1991, MNRAS, 252, 428 

El-Zant, A. A., Hoffman, Y., Primack, J., Combes, F.,& Shlos- 

man, I. 2004, ApJ, 607, L75 
El-Zant, A. A., Shlosman, I., & Hoffman, Y. 2001, ApJ, 560, 636 
Ettori, S., Fabian, A. C, Allen, S.W., & Johnstone, R.M., 2002, 

MNRAS, 331, 635 
Flores, R. A., & Primack, J. R. 1994, ApJ, 427, LI 
Fukushige, T., Kawai, A., Makino, J., 2004 ApJ 606, 625-634 
Gao, L., & White, S. D. M. 2007, MNRAS, 377, 5 
Gao L., Navarro J. F., Cole S., Frenk C. S., White S. D. M., 

Springcl V., Jenkins A., Neto A. F., 2008, MNRAS, 387, 536 
Gavazzi R., 2005, A&A, 443, 793G 

Gavazzi R., Fort B., Mellier Y., Pello R., Dantel-Fort M., 2003, 
A&A, 403, 11 

Gentile, G, Salucci, P., Klein, U., Vergani, D., & Kalberla, P. 

2004, MNRAS, 351, 903 
Gnedin, O. Y., Kravtsov, A. V., Klypin, A. A., & Nagai, D. 2004, 

ApJ, 616, 16 
Gottlober, S. & Yepes, G., 2007, ApJ 664:117-122 
Governato, F., et al., 2010, Nature 463, 203 

Graham, A.W., Merritt, D., Moore, B., Diemand J., and Terzic, 

B., 2006, AJ 132, 2701 
Graham, A.W., Merritt, D., Moore, B., Diemand, J., and Terzic, 

B., 2006, AJ 132, 2685 
Gregory S.A., Tiffit W.G., 1976, ApJ 205, 716 



Gustafsson, M., Fairbairn, M., & Sommer-Larsen, J. 2006, Phys. 

Rev. D, 74, 123-522 
Haarsma, D.B., et al. 2010, ApJ 713, 1037-1047 
Hayashi E. et al., 2004, MNRAS, 355, 794 
Hayashi, E., et al. 2004, MNRAS, 355, 794 
Hiotelis N., 2002, A&A, 383, 84 
Host, O., and Hansen, S.H., 2011, ApJ, 736, 52 
Hu, W., Kravtsov, A.V., 2003, ApJ 584, 702 
Hudson, M. J. & Ebeling, H. 1997, ApJ, 479, 621 
Hurley- Walker et al., 2011, arXiv: 1101.5912 
Icke V., 1973, A&A 27, 1 
Jing, Y. P., & Suto, Y. 2000, ApJ, 529, L69 
Kashlinsky, A., 1987, ApJ 312, 497 

Karachentsev I. D., Karachcntseva V. E., Huchtmeier W. K., 

Makarov D. I., 2004, AJ, 127, 203 
Karachentsev I. D., Kashibadze O. G, 2006, Astrophysics, 49, 3 
Keeton, C. R. 2001, ApJ, 561, 46 
Kelson et al. 2002, ApJ 576, 720 

Klypin, A., Kravtsov, A. V., Bullock, J. S., & Primack, J. R. 

2001, ApJ, 554, 903 
Klypin, A., Zhao, H-S., and Somerville R.S., 2002, ApJ 573, 597 
Kneib J. P., et al., 2003, ApJ, 598, 804 

Komatsu, E., Smith, K.M., Dunkley, J., Bennett, C.L., Gold, B. 

et al., 2011, ApJS, 192, 18 
Kowalski, M., Rubin, D., Aldering, G., Agostinho, R.J., Amadon, 

A. et al.,2008, ApJ, 686, 749 
Kravtsov, A. V., Klypin, A. A., Bullock, J. S., & Primack, J. R. 

1998, ApJ, 502, 48 
Lackner, C. N., Ostriker, J. P., 2010, ApJ 712, 88-100 
Laporte, C. F. P., White, S. D. M, Naab, T., Ruszkowski, M., 

Volker Springel, V., 2012. larXiv:1202.2357l 
Le Delliou M., Henriksen R. N., 2003, A&A, 408, 27 
Lewis, A. D., Buotc, D. A., Stocke, J. T., 2003, ApJ 586, 135L 
Limousin, M., ct al. 2008, A&A, 489, 23 
Lin C.C., Mestel L., Shu F.H., 1965, ApJ 142, 1431 
Lin, Y.T. & Mohr, J. J. 2004, ApJ, 617, 879 
Loeb, A., Peebles, P. J. E., 2003, ApJ 589, 29 
Lukic, Z., 2009, ApJ, 692, 217228 
Ma, CP., Chang, P. and Zhang, J.. larXiv:0907.3144l 
Mahdavi, A., Hoekstra, H., Babul, A., Sievers, J., Myers, S. T., 

& Henry, J. P. 2007, ApJ, 664, 162 
Mashchenko S., Couchman H. M. P., Wadsley J., 2006, Nat, 442, 

539 

McMillan, P. J., & Dchnen, W. 2005, MNRAS, 363, 1205 
Mellier Y., 1999, ARA&A, 37, 127 

Merritt, D., Navarro, J.F., Ludlow, A., and Jenkins, A., 2005, 
ApJ 624, L85 

Meneghetti, M., Bartelmann, M., Jenkins, A., Frenk, C, 2007, 

Volume 381, Issue 1, pages 171-186 
Miller C. J., et al., 2005, AJ, 130, 968 
Milosavljevic, M., & Merritt, D. 2001, ApJ, 563, 34 
Miralda-Escudc, J, 1995, ApJ 438, 514 

Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 
Moore, B. 1994, Nature, 370, 629 

Moore, B., Governato, F., Quinn, T., Stadel, J., & Lake, G. 1998, 
ApJ, 499, L5 

Morandi, A., Pedersen, K., & limousin, M., 2010, ApJ 713, 491 
Navarro ct al. 2010, MNRAS 402, 2134 
Navarro, J. F., et al. 2004, MNRAS, 349, 1039 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 
563 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 
493 

Newman, Andrew B.; Treu, Tommaso; Ellis, Richard S.; Sand, 

David J., 2011, ApJ 728, 39 
Newman, Andrew B.; Treu, Tommaso; Ellis, Richard S.; Sand, 

David J.; Richard, Johan; Marshall, Philip J.; Capak, Peter; 

Miyazaki, Satoshi, 2009, ApJ 706, 1078 



16 A. Del Popolo 



Nusser A., 2001, MNRAS, 325, 1397 

Oh, K. S. 1990, PhD thesis, Univ. California-Santa Cruz 

Oh S-H., Brook O, Govcrnato F., Brinks E., Mayor L., de Blok 

W. J. C, Brooks A., Walter P., 2010, AJ 142, 24 
Padmanabhan, T., 1993, Structure formation in the universe 

(Cambridge: Cambridge University Press) 
Peacock J. A., Heavens A.F., 1985, MNRAS 217, 805 
Peacock, J. A., & Heavens, A. F. 1990, MNRAS, 243, 133 
Peirani, S., Kay, S., & Silk, J. 2008, A&A, 479, 123 
Percival, W.J., Reid, B.A., Eisenstein, D.J., Bahcall, N.A., Bu- 

davari, T. et al., 2010, MNRAS, 401, 2148 
Pimbblet, K. A., MNRAS 411, 26372643 

Plionis, M., Tovmassian, H. M., Andernach, H., 2009, MNRAS 
395, 2 

Popesso, P., Biviano, A., Bohringer, H., & Romaniello, M. 2007, 
A&A, 464, 451 

Power, O, Navarro, J. F., Jenkins, A., Frenk, C. S., White, S. D. 
M., Springel, V., Stadel, J., & Quinn, T. 2003, MNRAS, 338, 
14 

Ricotti, M., and Wilkinson, M.I., 2004, MNRAS 353, 867 
Ricotti, M., 2004, MNRAS 2003, 344, 1237 
Ricotti, M., Pontzen, A., and Vicl, M., 2007, ApJ 663, 53 
Ryden, B. S., 1988, ApJ 329, 589 

Rix, H.-W., de Zeeuw, P. T., Cretton, N., van der Marel, R. P., 

& Carollo, C. M. 1997, ApJ, 488, 702 
Romano A., et al., 2010, A&A, 514, 88 

Romano-Diaz, E., Shlosman, I., Heller, C.,& Hoffman, Y. 2009, 
ApJ, 702, 1250 

Romano-Diaz, E., Shlosman, I., Hoffman, Y.,& Heller, C. 2008, 
ApJ, 685, L105 

Rood H.H., Page T.L., Kintner E.C., King I.R., 1972, ApJ 175, 
627 

Ryden, B. S., & Gunn, J. E. 1987, ApJ, 318, 15 
Salvador-Sole E., Solanes J.M., 1993, ApJ 417, 427 
Sand D. J., Treu T., Ellis R. S., 2002, ApJ, 574, L129 
Sand D. J., Treu T., Smith G. P., Ellis R. S., 2004, ApJ, 604, 88 
Sand, D. J., Treu, T., Ellis, R. S., Smith, G. P., Kneib, J.-P, 2008, 
ApJ 674, 711 

Schmidt, K.B., Hansen, S.H., and Maccio', A.V., 2008, ApJ 689, 
L33 

Schmidt, R. W., & Allen, S. W. 2007, MNRAS, 379, 209 
Schmidt, R. W., & Allen, S. W. 2007, MNRAS 379, 209 
Sharma, S., Steinmetz, M., 2005, ApJ 628, 21 
Schombert, J. M. 1988, ApJ, 328, 475 

Sikivie P., Tkachev 1. I., Wang Y., 1997, Phys. Rev. D, 56, 1863 
Simon J. D., Bolatto A. D., Lcroy A., Blitz L., 2003, ApJ, 596, 
957 

Simon J. D., Bolatto A. D., Leroy A., Blitz L., Gates E. L., 2005, 

ApJ, 621, 757 
Smith et al. 2001, ApJ 552, 493 

Smith, G. P., Kneib, J. P., Smail, 1., Mazzotta, P., Ebeling, H., 

& Czoske, O. 2005, MNRAS, 359, 417 
Smith, G. P., Kneib, J., Ebeling, H., Czoske, O., & Smail, f. 2001, 

ApJ, 552, 493 

Spano M., Marcelin M., Amram P., Carignan C, Epinat B., Her- 
nandez O., 2008, MNRAS, 383, 297 
Spekkens K., Giovanclli R., Haynes M. P., 2005, AJ, 129, 2119 
Spergel, D. N.; Verde, L.; Peiris, H. V.; Komatsu, E.; Nolta, M. 
R.; Bennett, C. L.; Halpern, M.; Hinshaw, G.; Jarosik, N.; 
Kogut, A.; 2003, ApJS 148, 175 
Stadel J., Potter D., Moore B., Diemand J., Madau P., Zemp M., 

Kuhlcn M., Quilis V., 2009, MNRAS, 398, 21 
Stadel, V. J., & Quinn, T. 2003, MNRAS, 338, 14 
Subramanian, K., Cen, R., and Ostriker, J. P., 2000, ApJ 538, 528 
Swatcrs, R. A., van Albada, T. S., van der Hulst, J. M., & Sancisi, 

R. 2002, A&A, 390, 829 
Swaters R. A., Madore B. F., van den Bosch F. C, Balcells M., 
2003a, ApJ, 583, 732 



Swaters R. A., Verheijen M. A.W., BcrshadyM. A., Andersen D. 

R., 2003b, ApJ, 587, L19 
Tonini, C, Lapi, A., and Salucci, P., 2006, ApJ 649, 591-598 
Treu, T., & Koopmans, L. V. E. 2002, ApJ, 575, 87 
Tyson, J. A., Kochanski, G . P., & dell' Antonio, 1. P., 1998, ApJ 

498, L107 

Umetsu, K., & Broadhurst, T. 2008, ApJ, 684, 177 
Umetsu, K., et al. 2009, ApJ, 694, 1643 

van den Bosch F. C, Robertson B. E., Dalcanton J. J., dc Blok 

W. J. G., 2000, AJ, 119, 1579 
Weinberg, M. D., & Katz, N. 2002, ApJ, 580, 627 
Whilcy, I. M., Aragon-Salamanca, A., Dc Lucia, G., et al. 2008, 

MNRAS, 387, 1253 
White, M., 2001, A&A Volume 367, 27-32 

Williams L. L. R., Babul A., Dalcanton J. J., 2004, ApJ, 604, 18 
Yang, X., Mo, H. J., & van den Bosch, F. C. 2008, ApJ, 676, 248 
Zappacosta, L., Buote, D. A., Gastaldcllo, F., Humphrey, P. J., 

Bullock, J., Brighenti, F., & Mathews, W. 2006, ApJ, 650, 

777 

Zitrin, A., et al. 2011, MNRAS DOl: 10.1111/j.l365- 
2966.2011.20155.x 



